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NON-EQUILIBRIUM THERMODYNAMICAL FRAMEWORK FOR 
RATE- AND STATE-DEPENDENT FRICTION 

P. VAN 1 ’ 2 , N. MITSUI 1 AND T. HATANO 3 


Abstract. Rate- and state-dependent friction law for velocity-step and heal¬ 
ing are analysed from a thermodynamic point of view. Assuming a logarithmic 
deviation from steady-state a unification of the classical Dieterich and Ruina 
models of rock friction is proposed. 


1. Introduction 


The rock experiments of sliding friction are understood by the so-called rate- 
and state-dependent friction laws. The equations of this laws unify the results 
obtained from two types of rock experiments; the first one is the time dependence 
of static coefficient of friction [I and the second one is slip velocity dependence of 
the dynamic coefficient of friction [2]. 

The properties of dynamic friction are the following 

(1) frictional coefficient in stable sliding conditions with a constant load-point 
velocity depends on the logarithm of the load-point velocity; 

(2) the magnitude of the instantaneous jump of the frictional coefficient de¬ 
pends on the change of the logarithm of the quotient of the corresponding 
load-point velocities; 

(3) the following evolution of the frictional coefficient to new value in stable 
sliding is also dependent on the instantaneous change of the load-point 
velocity; 

(4) oscillation occurs in some cases (e.g., large load-point velocity, polished 
surfaces, thin sand interface layer between the samples) (see e.g., H]). 

In healing experiments the properties of static friction are: 

(1) recovery magnitude is proportional to the logarithm of healing time; 

(2) larger velocity or larger elasticity increases the recovery magnitude of the 
static friction; 


These properties can be reproduced by using two classical equations. The first 
one is the constitutive law ([Tj) , expressing the relation between frictional coefficient 
/i and slip velocity V with an additional variable, called state variable 9. The 
second one is the evolution law (2]) expressing the time evolution of state variable 
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where r is the shear stress, er is the normal stress, /z* is the constant frictional 
coefficient for steady-state slip at reference slip velocity Vq, a and b are material 
parameters, L is the critical slip distance, and t is time. 

An important improvement is the evolution equation of Ruina [6]. 
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Experimental data of static friction is better reproduced by the equations of the 
Dieterich-law [5] (eqs. 0 and ©). and of dynamic friction by the equations of 
Ruina-law [6] (eqs. Q and ©)■ A comparison with experimental data is given in 
the works mm- In particular the Dieterich-law assumes time-relaxation and there¬ 
fore it is asymmetric for upward and downward jumps in displacement, contrary to 
the experiments. On the other hand the experimentally observed time dependent 
healing is properly reproduced by the Dieterich-law and does not reproduced by 
the Ruina-law. Thus another versions have been proposed (e. g ., [gam) in order 
to reproduce the experimental data better. However, none of them are completely 
satisfactory. 

Nakatani reformulated the Dieterich-Ruina law introducing a new variable m-- 

<4) s= i e>ti> (f)- 

Then the constitutive law is linearized with this variable: 

t fV\ 

(5) /z = -=/z*+alnl — 1+0, 


and the evolution laws of Dieterich and Ruina become 
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respectively. Nakatani suggested this modification together with a particular in¬ 
terpretation of the new variable as strength and interpreted the modified evolution 
equation of Dieterich in the framework of thermal activation theory. He did not 
investigate the modified form of Ruina law 0 . 


2. Thermodynamics of the frictional layer 

A thermodynamic approach of Mitsui and Van introduced a minimal model 
of rock friction by using permanent and recoverable parts of the total observed 
displacement as state variables in the spirit of continuum plasticity m- They 
interpreted the state variable of the rate- and state dependent friction law as the 
elastic part of the displacement. It was argued that the common origin of the 
constitutive and the evolution laws is originated in a thermodynamic framework. 
Their calculation of the entropy production was the following. 

The general setting is a sliding body on a horizontal surface with mass to. There 
are two forces that determine the motion of the body: the external force F ei and the 
damping force F d , due to friction (Fig. 0). The position of the body is denoted 
by x. The body is not considered completely rigid, however one assumes that one 
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Figure 1. Sliding thermomechanical body 

particular material point of the body characterize its instantaneous position. The 
equation of motion is 

(8) mx = F e - F d . 

Moreover, the work of the external force changes the energy of the body, E. 
Therefore 

(9) E = F e x. 

In this case thermodynamics requires that the damping force contributes only to 
the internal energy of the body. It is assumed that the external force accelerates the 
body and also that the body is deformable. In this homogeneous model, the defor¬ 
mation is expressed by the recoverable displacement, r. Therefore the kinetic and 
and elastic energies of the body are distinguished. This particular interpretation 
from mi is not necessary, r may denote a general internal variable. 

The internal energy, U, is the difference of the total energy, E, the kinetic energy 
and a quadratic contribution of the internal variable. This form follows from the 
condition of thermodynamic stability in the state space El- 

( 10 ) U = E -m^- k r Y , 

where k is a parameter like in elasticity. One assumes a particular kinematic condi¬ 
tion in order to introduce an interpretation of r. When the instantaneous position 
of the body is the sum of a permanent and a recoverable displacements then a con¬ 
venient method of their distinction is an additive separation of the displacement 
rates: 

( 11 ) 


V = x = r + z 
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where V is the rate of the position x, and z is the rate of the permanent dis¬ 
placement. This rate type kinematical condition is convenient when distinguishing 
permanent and recoverable changes. © is analogous to the condition used in plas¬ 
ticity for the distinction of plastic and elastic strains (see e.g., |l.'ll I14jl. However, 
in friction the internal variable is not necessarily identical to the recoverable strain. 

The dissipation can be calculated by the entropy balance, assuming that the 
entropy is the function of the internal energy only: 

(12) S{U) = ^ (F e x - mVV - fcrr) >0 => TS = F d V - krr > 0. 


The damping force and also the rate of the internal variable f are the constitu¬ 
tive quantities to be determined in accordance with the requirement of nonnegative 
entropy production. We connect the first term to friction interaction and the sec¬ 
ond term to frictional healing. In the vicinity of thermodynamic equilibrium, where 
both forces and fluxes are zero, the usual linear relationship is a consequence of La¬ 
grange’s mean value theorem [1511X5] , V and r are thermodynamic state variables, 
therefore they can be considered as thermodynamic forces, while /r and nr are to 
be determined constitutively, they are thermodynamic fluxes. The standard linear 
approximation results in equations that reflect well the thermodynamic admissibil¬ 
ity of both velocity weakening and strengthening, however, did not incorporate the 
observed direct effect (logarithmic relaxation) [TTI . 

However, in velocity-step experiments dynamic friction is a steady state phe¬ 
nomenon, we are looking for a deviation of the frictional coefficient from a fixed 
value, fjj o, at a given Vo reference velocity. On the other hand the internal variable 
is expected to be zero, when its time derivative is zero, in thermodynamic equi¬ 
librium. Therefore we face to a mixed, partial steady state, partial equilibrium 


situation. However, the dissipation inequality (121 may not characterize the devi¬ 
ation from steady state. What we need is an estimation of the deviation of the 
entropy production. 

In our framework the steady state means a constant displacement rate, denoted 
by Vb- Therefore the internal energy should be determined accordingly 


(V-V 0 y . r 


2 r 2 

— k- 


(13) U = E 0 - to- 2 .. 2 

where E 0 is the energy of the body moving with the veolcity V)f] Then the calcu¬ 
lation of the entropy production results in 


(14) TS(U) = F e V - m(V - V 0 )V - krr = F e V 0 + F d (V - V 0 ) - krr > 0. 


We can introduce a friction specific entropy production with Amonton’s law: 

(15) S = S = hqVq + A/i(V - V 0 ) - Krr > 0, 

■F n-A. 

where = F e /(AF n ) is the external shear stress divided by the normal force, 
A fi = Fd/(AF n ) is the frictional stress divided by the normal force and k = 


^The proper relation of kinetic and internal energies in a Galilei relativistic framework is a 
delicate question from a thermodynamic point of view. A detailed treatment of objective, frame 
indpendent thermodynamic modeling in case of single component fluids is given in m- A similar 
approach is required in our case, too 
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Now, it is straightforward to introduce a linear approximation for the increment 
of entropy production, as for the near equilibrium situation. However, in the fol¬ 
lowing we apply a different starting point. Our fundamental assumption is that the 
leading term of the deviation from the steady state is logarithmic , while it is linear 
around equilibrium. We can formulate this hypothesis analogously to the classical 
exploitation of the entropy principle introducing the concept of incremental entropy 
production. We require that it is minimal at the steady state 

(16) AE = E - T, st eady = E - n 0 V 0 = A/iln ~ ^ > 0 

We will call this hypothesis the principle of minimal incremental entropy production. 

The required minimality ensures the asymptotic stability of the steady state if 
the evolution equations are constructed accordingly. In this respect the principle is 
similar to the role of the second law near to the equilibrium. On the other hand this 
hypothesis is a modification of the requirement of nonnegative excess entropy pro¬ 
duction of Glansdorff and Prigogine [IS]. The difference of the Prigogine-Glansdorff 
requirement is the logarithmic deviation instead of a linear one. (The possibility of 
partial steady state, a specific property of friction, is important, too). Along with 
them we want to emphasize that the inequality here is not a law of nature, it is 
regarded as a convenient stability assumption mm- 

An example of similar logarithmic deviation is the thermodynamics of chemical 
kinetics, where the entropy production is a product of the chemical affinity and the 
reaction rate, but the Guldberg-Waage kinetic equations introduce an exponential 
relation m- It is remarkable that chemical equilibrium is considered as a steady 
state from a thermodynamic point of view, when forward and backward reactions 
are properly distinguished [22] 123] . Our formula corresponds a simple monomolec- 
ular reaction, which is actually equivalent to a single internal degree of freedom 
|21| . In friction the logarithmic deviation may be further motivated according to 
this chemical analogy [ID] . 

Finally we remark, that the logarithmic form can be directly derived assum¬ 
ing that the internal energy is modified by a logarithmic velocity dependent term 
instead of a quadratic one U = E 0 — mV(Log(V/Vo ) — 1). 


(161 is similar to the usual entropy production in many respect. First of all the 


two terms are zero in the reference steady state and we assume that A/r is a consti¬ 
tutive function of the logarithmic deviation, a function of the thermodynamic state 
variable. Therefore, in case of smooth functions the Lagrange mean value theorem 
ensures a linear homogeneous relationship between these constitutive quantities and 
the related thermodynamic forces also in this mixed steady-state equilibrium case. 


The linearization of (16) result in the following expression: 


(17) 

(18) 


A/.t = l\ In 
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We will call © -( |18| ) as thermodynamical aging law. © is identical to the 
Nakatani form of the constitutive law ([5j) and (18) is similar to |7|. 

The coefficient matrix may depend on the thermodynamic forces, in particular 
h, h, hi and Z 2 i may depend on the state variables V and r [Mj. In the following 
we assume a strict linear relationship, when the coefficient matrix is constant. It is 
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remarkable, that there are no reasons to assume symmetry or antisymmetry of the 
matrix. The conditions of Onsagerian statistical background cannot be introduced 
without a particular interpretation of the internal variable (more detailed arguments 
are given in [251 W. ■ 


3. Different mechanisms of different relaxations 


According to the experimental observations in case of velocity step and healing 
experiments, sometimes the internal variable changes when the surfaces slip, and 
sometimes its evolution is seemingly independent of the relative motion of the 
surfaces. This distinction is connected to the detailed mechanism of friction, and 
requires an extension of the modeling framework. In the following we will investigate 
the question of slip related internal variable evolution. 

Ruina [ft applied this assumption directly to the relaxation. In our case that 
requires the modification of (18), assuming that the slip is what makes change of 
the internal variable and the rate of the variable is a consequence. In this case 


the time derivative in (18) is substituted by the space derivative and the rate is 


obtained as a consequence: 
(19) 


dr 

dt 


dr 

dx 


r 

V' 


Performing this substitution in (18) leads to 


( 20 ) 


r = l 21 V ln( ^ ) -VhKr. 


This equation, together with will be called thermodynamical slip law. One can 
see, that the Nakatani transformed Ruina-law <m>, 0 can be obtained if l\ = a, 
l \2 = — 1/k, l 2 \ = —h/L and l 2 = 1 /{kL). The number of phenomenological 
parameters is increased by one, from three to four, compared to the original Ruina 
theory. It is because k appears only as a multiplier of the cross coefficients l\ 2 and 
l 2 i. The slip governed modification represents a particular quasilinear form of the 
strictly linear relations of the thermodynamic aging law (fl7[)-(|l8t 


The slip condition of Ruina, expressed by (19) assumes that V is the relative 


surface velocity, related to the permanent part of the displacement. When intro¬ 
ducing a distinction between a permanent and a recoverable parts of the apparent 
displacement beyond the difference of the load point an relative surface velocities, 
one may expect, that only the permanent part contributes to the entropy produc¬ 
tion, to the evolution of the internal variable. For example an interpretation of the 


2 In our thermodynamical framework a slip related change may lead to further consequences. 
In order to keep the integrity of the thermodynamic considerations, the calculation of the entropy 
rate may be substituted by the calculation of the slip related entropy change. 

For example, with the internal energy, © the derivative of the entropy by the displacement 
will be the following: 


( 21 ) 


dS 

dx 


1 

T 



+ F d 


V — Vo 
V 



Then we proceed assuming logarithmic increment and obtain © and ( |20| ) as a consequence. 
However, slip and displacement are not the same, the calculation of slip related changes should 
distinguish between the permanent and recoverable parts. Here we do not analyse this possibility 
further, we accept the approach of Ruina at this point. 
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internal variable as recoverable displacement with the condition (11) leads to the 
following permanent displacement: 


X per = / z(s)ds = X — T 
Jt 0 

Particular mechanisms require that the internal variable influences the slip. For 
example when the internal variable is connected to deformation of surface irreg¬ 
ularities, then this conclusion is straightforward. Therefore we assume in gen¬ 
eral, that the internal variable directly influences the displacement and the reduced 
part is what influences the (incremental) entropy production. Hence we introduce 
x re d = x — ar, where 0 < a < 1 is the factor of slip reduction. If a = 1, then 
the internal variable can be interpreted as the recoverable part of the displacement 

HU 

Therefore the evolution equation of the internal variable is obtained by substi¬ 
tuting the time derivative with the slip related change of the internal variable, as 
follows: 

dr dr f 

dt d(x — ar) V — ar 

For the sake of simplicity a is constant. Then a simple rearrangement leads to the 
following evolution equation of the internal variable: 


( 22 ) 


r =V 


hi 
1 + a 


In (^) - h^r 
[hi ln(^) -hnr} 


Together with © we will call this equation thermodynamical friction law. This 
is an other particular quasilinear form of the thermodynamic aging law ©-@. 

In the above thermodynamical models there is a direct effect and the conventional 
step test parameters are a = 1 1 and b = lv \ l21 . In the following we compare the 
performance of the obtained thermodynamic models with the classical models and 
experiments. 


4. Velocity step tests 

A comparison of the different rate- and state dependent friction laws is shown 
in Figure [2] The velocity weakening experiment is modeled with the following 
parameters: 

po = 0.6, Vo = 1pm/s, Vi = 10 pm/s, L = 20 pm, a = 0.015, and b = 0.02. 

Dieterich and Ruina models shown by a solid thin lines. The parameters of the 
thermodynamic friction model are calibrated to give the proper step conditions 
and relaxation speed: l\ = a, h = 1/L, l\2 = b, hi = h- There are two additional 
parameters, the static recovery strength n, analogous to a spring constant, and the 
factor of slip reduction a. The dotted curve runs exactly over the line obtained by 
the Ruina model, because thermodynamic slip model recovers the Ruina case when 
k = 1 and a = 0. One obtains highly asymmetric relaxation curves choosing higher 
values of the slip reduction parameter a and with a proper choice one calculates 
curves that run close are pretty close to the Dieterich model. E.g. the thick 
dashed curve was calculated with n = 0.85 and a = 8. Moreover, one can obtain 
symmetric relaxation curves close to either the down or up relaxation curves of the 
Dieterich model with an appropriate choice of the parameters. For example the 
thick dotdashed line was calculated using k = 0.7 and a = 2. 
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The thermodynamic aging law produces asymmetric curves, similar to the Di- 
eterich law. 





s[//m] 


Figure 2. Simulation of a velocity step test with different friction 
laws. The Dieterich and Ruina laws are drawn by solid thin lines. 
The thermodynamic friction law leads to the dotted line with k = 1 
and a = 0 parameters running over the line of Ruinalaw. The 
dashed line with parameters k = 0.85 and a = 8 runs close to the 
Dieterich law. For the dotdashed line the parameters are k = 0.7 
and a = 2. 


5. Healing 

The interpretation of healing is contradictory, therefore we have chosen the ex¬ 
periments and strategy of Beeler an Tullis for demonstration. They per¬ 

formed healing experiments with different machine rigidity. One of them k = 0.002 
was the natural elasticity of the experimental device, in this case the load point 
velocity was V 0 = 1/rni/s. With the k = 0.074 servocontrolled value the load point 
velocity was Vq = 0.316 fim/s. In both cases the rock parameters were L = 3 /im, 
a = 0.009, and b = 0.008. The data points in the first and second cases are shown 
on figure [3] by circle and rectangles, respectively. The simulation used the exper¬ 
imental rock parameters of Beeler et al. (271 and the additional parameters were 
chosen as a = 0.75 and n = 0.9. Then we have obtained the dashed curve for the 
k = 0.074 case and the dotted one for k = 0.002. 

Figures [4] and [5lshow the effect of changing the initial velocity and the machine 
rigidity. In figureH] the simulation of the healing experiment with the more rigid 
machine is shown, where the parameter values are k = 0.002, Vq = l/x?n/s, L = 
3 /im, a = 0.009, b = 0.008, a = 0.75 and n = 0.9 (dotted curve in figure [3|. 
Increasing the velocity to Vq = 2 [im/s pushes the curve upward and parallel to 
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Figure 3. Simulation of the healing experiments of Beeler et al 
m- The experimental data for the k = 0.002 case is shown by 
the big circles and for the k = 0.074 case by the rectangles. The 
dotted curve is calculated by the thermodynamic friction model 
for k = 0.002. The k = 0.074 parameter resulted in the dashed 
curve. 


the original one, shown with the dashed curve. In figure [5] the effect of softening 
is demonstrated, the dashed curve is calculated with k = 0.02. The dotted curve 
is identical to the one in figure |4j The increase of the healing effect qualitatively 
corresponds to the experimental observations. 

6. Summary and conclusion 

A non-equilibrium thermodynamic model with a single internal variable of rate- 
and state-dependent friction was proposed. The model introduced the following 
basic assumptions: 

- The deviation from the steady state is logarithmic. 

- The change of internal variable is due to a reduced slip, and the reduction 
is proportional to the value of the internal variable. 

The obtained thermodynamic friction model generalizes the well known Dieterich 
and Ruina laws with two additional parameters. 

- For velocity step tests it interpolates between the Dieterich and Ruina laws. 
The form of the relaxation depends on the additional parameters. One can 
obtain symmetric up and down relaxation curves that are either close to 
the curves of Ruina relaxation or are close to the Dieterich type relaxation, 
but with less apparent linear part. 

Simulations show promising results for healing experiments, both quanti- 
tavely and qualitatively. 
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Figure 4. The effect of larger steade state velocity for healing. 
Vq = 1 /j,m/s, dotted curve, Vq = 2 /j.m/s, dashed one. 
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Figure 5. The effect of softened machine for healing, k = 0.002, 
dotted curve, k = 0.02, dashed one. 
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